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Abstract 

We have studied the superconducting, magnetic and charge correlation func- 
tions and the spin excitation spectrum in the doped two chain Hubbard model 
by projector monte carlo and Lanczos diagonalization methods. The exponent 
of the interchain singlet superconducting correlation function, 7 , is found to 
be close to 2.0 as long as two distinct non- inter acting bands cross the Fermi 
level. Magnetic and charge correlation functions decay more rapidly than or 
as fast as the interchain singlet superconducting correlation function along 
the chains. The superconducting correlation in the doped two chain Hubbard 
model is the most long range correlation studied here. Implications of the re- 
sults for the possible universality class of the doped two chain Hubbard model 
are discussed. 

PACS numbers: 71.27.-ha,74.20.Mn.,75.40.Mg,75.40.Gb 
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I. INTRODUCTION 



In relation to possible superconductivity in the two dimensional Hubbard model 
the two chain Hubbard model and/or the two chain t-J model, which may be realized 
in oxygen deficient S'r„_iC'U„+i02„ have attracted much interest recently. Although 
superconductivity has not yet been observed in this material, it may be worthwhile to study 
these models because they may make it possible for us to understand dimensional cross- 
over between one and two dimensions. Information obtained in them may be of help in 
understanding the two dimensional models. The two dimensional models are more difficult 
to study than the one dimensional models and the two chain models. Though direct studies 
of the two dimensional models are also very important, here we follow the argument in 
Refs. and study the doped two chain Hubbard model, which is defined by the following 
Hamiltonian: 



H = -h E (4c,. + H.C.) + UY: nfr^nfri 

-t^Y.i4A. + H.c.). (1) 

We take ty = 1.0 and adopt it as the unit of energy hereafter. Diagonalization studies 
of the doped two chain models suggest that the superconducting correlation function may 
be rather long range. [|7|,p!0| The cluster size accessible by the Lanczos method is obviously 



too small compared with the superconducting correlation length in these models to make a 
convincing conclusion. The problem should be examined with more sophisticated methods 
in computational physics without any approximation and with sufficiently large system size. 

The recently developed density-matrix real space numerical renormalization group 
method has been applied to study the long range superconducting correlation in the doped 
two chain Hubbard model. The authors concluded that there is an enhancement of the 
superconducting correlation function over the non-interacting case on up to as large as 32 x 2 
system sizes when tj_ = 1.5, U = 8.0 at N^/Ns = 0.875 filling, where Ng and Ng are the num- 



ber of electrons and the number of sites, respectively. They also have found persistence of 
the spin gap after doping with the same parameter values. 

Though their method seems to be quite useful at least in some of one dimensional spin 
systems, we take here a different approach to study some interesting phases of the doped two 
chain Hubbard model. For most of this article we adopt the projector quantum monte carlo 
method. While the method is quite useful to our model in the small-U region, the method 
currently has difficulty in the large-U region due to the fermion sign problem. Our strategy 
here is to derive reliable results in the small-U region, which still has not been explored as well 
as to get some preliminary information in the large-U region. The similarity or dissimilarity 
of the two regions should be of great interest. So long as each numerical method has both 
pros and cons, the author believes that the problem should be approached with various 
numerical techniques. 

II. STATIC AND DYNAMIC CORRELATION FUNCTIONS 

Here we study singlet superconducting < O^O] >, magnetic < S?Sj > and charge 

< riirij > — < n >ave correlation functions and dynamic spin correlation function S^^(k,Oi)) 
in the ground state of the doped two chain Hubbard model. We have employed the projector 
monte carlo method to study the interchain singlet superconducting correlation function 

< OjO] >, where 

0]-^(c3t4-44)' (2) 
magnetic < SfSj > and charge < riirij > — < n >^^g correlation functions, where 

St ^ ^{rii^ - nil) (3) 

and 

rii = riii + riii (4) 

and the Lanczos and the recursion methods to study S^^{k,u!): 

3 



S''''{k,uj) = -——ImG{k,Eo + uj + ir]), (5a) 
Nn 

Gik,x) = {iPoliSiyix-Hy'Sll^o). (5b) 

We have adopted periodic boundary condition throughout this article. 36 x 2 and 24 x 2 
lattices were used for projector monte carlo calculations. A 6 x 2 lattice was used for Lanczos 
and recursion calculations. The projecting time r = 10.0 and r = 1.5 in units of ty were 
used for U = 2.0 and U = 6.0, respectively, in the projector monte carlo calculations. We 
have adopted the unrestricted Hartree-Fock (UHF) wavefunction as a trial function of the 
projector monte carlo calculation here. Details of these methods should be looked up in the 
previous papers and references cited therein. The two chain Hubbard model with various 
parameter sets at 0.833 filling {N^/Ns = 0.833) is studied here. The singlet superconducting 
correlation function in the non-interacting case is given by: 

< 0,0] >u=o= - cos{2k0)l) - cos{2kf{7c)l)], (6) 

where fc/(0) = cos~^(t± + /i)/2, kfi^n) = cos~^(t± — /i)/2 and 7 = 2, when U = 0.0. The 
upper bound of < OjOj > in the non-interacting case is given by 0.05-R~-^, where / and Rij 
are distance between the i-th and j-th sites. 

First of all, we study long-range behaviors of some of static correlation functions in 
small-U region of the model. < OiO] > when = 1.5, f/ = 2.0 and = 2.5, U = 2.0 
were calculated and they are plotted along with the upper bound of the non-interacting 
case in Fig. |l[ It should be noted that while two distinct bands cross the Fermi level when 
t± = 0.5, U = 0.0 and t_i = 1.5, U = 0.0, one of them is raised in energy and does not cross 
the Fermi level when t± = 2.5, U = 0.0. The correlation function with t± = 2.5, U = 2.0 
is strongly suppressed from the non-interacting one, which does not contradict the result 
obtained from the more general superconducting correlation function. |^ The suppression 
may originate from the separation of the bands. The correlation function with t± = 1.5, U = 
2.0 is somewhat enhanced or remains comparable to the non-interacting correlation function. 
This parameter region should be examined more carefully. 
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The superconducting correlation function < OjO] > with t± = 1.5, f/ = 2.0 was re- 
plotted along with the upper bound of the non-interacting case on a log-log scale in Fig. 

H. The correlation function has inversion symmetry at the inversion center: |i — j| = 17.5. 
We show only the left hand side of the function. Though the correlation function itself is 
somewhat enhanced over the non-interacting one, the exponent of the correlation function 
7 remains close to that of the non-interacting case: 2.0. 

At this stage, we compare the long range behaviors of the interchain singlet supercon- 
ducting, magnetic and charge correlation functions. We have calculated these correlation 
functions with t± = 1.5, U = 2.0 in 24 x 2 lattice and plotted them in Fig. ^. As mentioned 
earlier, < OjOj > scales as 0.05i?^-^. The data points of the magnetic < S^Sj > and charge 
< riiTij > — < n >^^g correlation functions are found between ±0.05-R~-^ and that means 
the magnetic and the charge correlation functions decay more rapidly than or as fast as the 
interchain singlet superconducting correlation function. The interchain singlet supercon- 
ducting correlation function is the most long range correlation function among these three 
correlation functions. This is quite different from finite-U one dimensional Hubbard model. 

The superconducting correlation function < OjOj > with tx_ = 1.5, f/ = 6.0 and t±_ = 
0.5, U = 2.0 were plotted along with the upper bound of the non-interacting case in Fig. 
^. < OiOj > for both the two parameter sets scale as 0.05-R~-^ like < OjO] > for t± = 

I. 5, U = 2.0 does. The oscillatory part for t± = 1.5,U = 6.0 is more enhanced than that for 
t± = 1.5, U = 2.0. The oscillatory part with t± = 0.5, U = 2.0 is largely reduced compared 
with that of t_i_ = 1.5, U = 2.0. The fermion sign problem in our present version of projector 
monte carlo method does not allow us to use larger projecting time than r = 1.5 when t± = 
1.5, U = 6.0 at 0.833 filling. The result with U = 6.0 should be taken preliminary. Whether 
or not the large-U region and small-U region can be classified into the same universality 
class (possibly with a spin gap) or not may be a subject of future interest. 

We are not able to decide if the spin gap persists after doping or not from our data on 
the correlation function because of the size of the error bars and the oscillatory behavior. 
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The authors of Refs. |17|,|9|JTT| beheve that the spin gap persists based on their numerical 



results. If the spin gap persists after doping, one of the possibilities is that our model 
is classified into the Luther-Emery (LE) liquid phase with exponential decay of the 2k p 



magnetic correlation function, (see Table II of Ref. |]T2[) As the exponent of power-law 
decay of < OjO] > , 7 , is estimated to be about 2, the exponent of 2kF charge correlation 
function should be about 0.5 in the LE liquid phase. This leads to a conclusion that the 
2kF charge correlation is the most long-range. Clearly, our data with U = 2.0 (Fig. is 
not consistent with this LE picture. Some other exotic possibilities consistent with the spin 
gap should be pursued. [|lT| In spite of slow convergence of lnb{2kF, L/2) versus InL in the 



Fig. 5 of the previous communication [Q, we may not be able to exclude a possibility of the 
Tomonaga-Luttinger (TL) liquid phase in the small-U region, however. 

The dynamic spin correlation function was calculated. The elementary spin excitation 
spectrum was obtained by tracing out u , giving the first peak of S^^{k, cu), where k = {k\\, k±) 
with defined along the chains and k±_ defined vertical to the chains. It should be noted that 
the spectrum can be also obtained from the energy eigenvalues in subspaces with definite k 
values. Both the k± = and k± = it branches of the spectra are plotted in Figs. |^. Three 
different parameter sets; = 0.5, U = 2.0, t± = 1.5, U = 2.0 and t± = 1.5, U = 8.0 have 
been used. The spectra at half-filling {Ne/Ng = 1.0) are shown in the insets. Spectra of 
the three parameter sets at half-filling are isomorphic to each other irrespective of values 
of the parameters in both the k± = and the k^ = ir channels. In the k^ = channel, 
k± = 7T is no longer one of the lowest energy states unlike in the one dimensional half-filled 
Hubbard model. At 0.833 filling, spectra for t±_ = 1.5, U = 2.0 and t± = 1.5, U = 8.0 are 
isomorphic to each other in both the k± = and the k± = n channels. The spectrum for 
t± = 0.5, U = 2.0 is quite different from them. The spin dynamics in the doped two chain 
Hubbard model changes with t^. The change occurs between t±_ = 0.5 and t±_ = 1.5. We 
do not observe the change in the half-filled state. The reduction of the oscillatory part of 
< OiO] > for ti_ = 0.5, U = 2.0 may originate from the change in spin dynamics. The change 
has no effect on 7. In the non-interacting case, we have two separated bands when t± 7^ 0. 



Both of the two bands are partially filled in both the two cases studied here; t± = 0.5 and 
t± = 1.5. It seems unlikely that the change in the spin dynamics is related with the number 
of partially filled bands. The Fermi wave vector of the second band (higher in energy) kp2 in 
the noninteracting case is 7r/3, which is closest to 7r/2 in this system size. The spin dynamics 
when t± = 0.5 and U = 2.0 at 0.833 filling may be related with the phase with two gapless 
spin modes and one gapless charge mode found by using weak coupling renormalization 
group techniques (the C1S2 phase in the notation of the authors of Ref. [Il3|). 

It should be noted that the smallest energy differences with the ground state are given at 
(fc||, k±) = (0.177r, 0.0), (0.507r, 0.0), (O.GTtt, 0.0) and (0.837r, 0.0) for = 0.5, U = 2.0 and at 
(O.ITtt, vr), (0.507r, vr), (O.GTvr, vr) and (0.837r, vr) for = 1.5, U = 2.0 and = 1.5, U = 8.0, 
respectively. The smallest energy differences for these parameter sets are 0.08, 0.10 and 0.37, 
respectively. The possible spin gap for U = 2.0 is less clear than that with U = 8.0 in this 
small system size. 

III. CONCLUSIONS 

We have calculated the interchain singlet superconducting, magnetic and charge corre- 
lation functions and the spin excitation spectrum in the doped two chain Hubbard model. 
The exponent of the interchain singlet superconducting correlation function, 7 , is found 
to be close to 2.0 as long as two distinct non-interacting bands cross the Fermi level. Our 
data is not consistent with the LE picture. Some other exotic possibilities consistent with 
the spin gap should be pursued. Of the three correlation functions, the interchain singlet 
superconducting, magnetic and charge correlation functions, the interchain singlet super- 
conducting correlation function is the most long range correlation function in the doped two 
chain Hubbard model. This is quite different from finite-U one dimensional Hubbard model. 
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FIGURES 

FIG. 1. The superconducting correlation function < OiOj > at 0.833 filling on a 36 x 2 lattice 
is plotted versus \i — j\- Circles and squares are the data obtained with t± = 1.5, U = 2.0 and 
t± = 2.5, U = 2.0, respectively. The upper bound of < OjOj >?7=o: O.ObRj^j is plotted as well, 
where Rij is distance between sites in the rung. 

FIG. 2. The superconducting correlation function < OiOj > at 0.833 filling on a 36 x 2 lattice 
is plotted versus |^ — j| on a log-log scale. Closed circles connected with real line are the data 
obtained with t± = 1.5, U = 2.0. A dotted line is upper bound of < OiOj >u=o' 0.05\i — 
The correlation function has inversion symmetry at the inversion center: \i — j\ = 17.5. We show 
only the left half of the function. 

FIG. 3. The superconducting < OjOj >, magnetic < SfSj > and charge < riirij > — <n >ave 
correlation functions at 0.833 filling on a 24 x 2 lattice with t± = 1.5, U = 2.0 are plotted versus 
\i — j\. Circles, squares and crosses are the superconducting, magnetic and charge correlation 
functions, respectively. The functions ±0.05i?~^ are also plotted as guides. 

FIG. 4. The superconducting correlation function < OjOj > at 0.833 filling on a 36 x 2 lattice 
is plotted versus \i — j\. Circles and squares are the data obtained with t± = 1.5, U = 6.0 and 

t —9 

t± = 0.5, U = 2.0, respectively. The upper bound of < OjOj >u=Q' 0.05i?jj- is plotted as well, 
where Rij is the distance between sites in the rung. 

FIG. 5. The spin excitation spectrum of a 6 x 2 two chain Hubbard model at 0.833 filling. 
The k± = and k± = ir channel are plotted in (a) and (b), respectively. Open squares, closed 
squares and closed circles are the spectra obtained with t± = 0.5, U = 2.0, t± = 1.5, U = 2.0, 
and t± = 1.5, U = 8.0, respectively. The horizontal axis is the wavenumber in units of tt. The 
corresponding spectra at half-filling {Ne/Ng = 1.0) are depicted in the insets. 
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